Mixing-demixing in a trapped degenerate fermion-fermion mixture 
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We use a time-dependent dynamical mean-field-hydrodynamic model to study mixing-demixing 
in a degenerate fermion-fermion mixture (DFFM). It is demonstrated that with the increase of 
interspecies repulsion and/or trapping frequencies, a mixed state of DFFM could turn into a fully 
demixed state in both three-dimensional spherically-symmetric as well as quasi-one-dimensional 
configurations. Such a demixed state of a DFFM could be experimentally realized by varying an 
external magnetic field near a fermion-fermion Feshbach resonance, which will result in an increase 
of interspecies fermion-fermion repulsion, and/or by increasing the external trap frequencies. 

PACS numbers: 03.75.Ss 
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I. INTRODUCTION 

Recent successful observation of degenerate boson- 
fermion mixture (DBFM) of trapped alkali-metal atoms 
by different experimental groups 0, 0> 0, 0] initiated in- 
tensive experimental studies of different novel phenom- 
ena 0, 0, 0- Among these experiments, apart from 
the study of DBFM 6 - 7 Li Q, 23 Na- 6 Li and 87 Rb- 40 K 
[11113, there have been studies of degenerate fermion- 
fermion mixtures (DFFM) of two components of 40 K 1] 
and 6 Li 0, 0. The collapse in a DBFM of 87 Rb- 40 K 
atoms has been observed [E[ and studied by Modugno 
et al. and by the present author ^(j and more recently by 
Ospelkaus et al. 8] . There have also been investigations 
on the formation of bright and dark jl2j solitons in 
a DBFM. 

Several theoretical investigations 0, 0, of a 
trapped DBFM studied the phenomenon of mixing- 
demixing when the boson-fermion repulsion is increased. 
For a weak boson-fermion repulsion at K both the Bose- 
Einstein condensate (BEC) and the degenerate Fermi gas 
(DFG) have maxima of probability density at the cen- 
ter of the harmonic trap. However, with the increase of 
boson-fermion repulsion, the maximum of the probabil- 
ity density of the DFG could be slowly expelled from the 
central region. With further increase of boson-fermion re- 
pulsion, the DFG could be completely expelled from the 
central region which will house only the BEC. This phe- 
nomenon has been termed mixing-demixing in a DBFM. 
The phenomenon of demixing has drawn some attention 
lately as in a demixed state an exotic configuration of the 
mixture is formed, where there is practically no overlap 
between the two components and one can be observed 
and studied independent of the other. It has been ar- 
gued [l5| that such a demixed state in DBFM should 
be possible experimentally by increasing the interspecies 
scattering length near a Feshbach resonance. In view of 
this it is of interest to see if, in addition to a DBFM, such 
a demixed state appears in a DFFM. 
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The purpose of this paper is to study and illustrate this 
mixing-demixing phenomenon in a trapped DFFM us- 
ing a coupled time-dependent mean-field-hydrodynamic 
model inspired by the success of a similar model in the 
investigation of fermionic collapse ^} an( l bright 
and dark solitons in a DBFM. The conclusions of 
the study on bright soliton are in agreement with a 
microscopic study [lfj], and those on collapse jlCj are in 
agreement with experiments 0, @ ■ There have been prior 
suggestions of mixing-demixing in a trapped DFFM upon 
an increase of interspecies repulsion. Roth and Feldmeier 
[l7| discuss the possibility of demixing mathematically 
using a, energy density functional approach. Amoruso et 
al. [lSj demonstrate demixing numerically using model 
hydrodynamic flow equations. In contrast to these pre- 
vious time-independent studies for stationary states, the 
present study relies on a time-dependent formulation and 
demonstrates demixing also for an increase of the trap- 
ping frequencies and should be considered complimentary 
to these previous investigations. The present suggestion 
of demixing with an increase of trapping frequency is eas- 
ier to implement experimentally than increasing the in- 
terspecies scattering length using a Feshbach resonance. 

We look for mixing-demixing in a trapped DFFM 
in the spherically-symmetric as well as quasi-one- 
dimensional configurations. For a weak fermion-fermion 
repulsion, we find the formation of a fully mixed state 
with both the fermion clouds occupying the central re- 
gion of the trap. If the fermion numbers of the two types 
of fermions are largely different, a fully demixed state 
of DFFM appears with the increase of fermion-fermion 
repulsion, when the component with smaller number of 
fermions is completely expelled from the central region 
of the trap which is occupied only by the component 
with the larger number of fermions. This demixing also 
appears for increasing trap frequencies which increases 
the fermion density and hence fermion-fermion repulsion 
thus resulting in a demixing. To the best of our knowl- 
edge this is the first study of mixing-demixing in a DFFM 
using a time-dependent formulation capable of investigat- 
ing noncquilibrium states. Previous studies of demixing 
|l7l 118] employed time-independent formulations appro- 
priate for stationary states only. 

In Sec. II we present the mean-field-hydrodynamic 



2 



model we use in our investigation in the spherically- 
symmetric as well as quasi-one-dimensional configura- 
tion. The numerical results are presented in Sec. Ill 
where we illustrate demixing with the increase of inter- 
species repulsion and of the harmonic trapping potential. 
Finally, in Sec. Ill a brief summary of our findings are 
given. 



II. MEAN-FIELD-HYDRODYNAMIC MODEL 

A proper description of a DFFM should employ (lHIHl ] 
a fully antisymmetrized microscopic Slater determinant 
wave function for each component. However, a simpli- 
fied mean-field-hydrodynamic Lagrangian for a DFG has 
been used successfully to study a DBFM 0, , which 
we shall use in the present investigation. The virtue of 
the mean-field model over a microscopic description is 
its simplicity and good predictive power. To develop a 
set of practical time-dependent mean-field-hydrodynamic 
equations for a DFFM, we consider the following La- 
grangian density [n], [T^] 
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where rrij is the mass of component j{— 1,2), a 
complex probability amplitude, nj — \^j\ 2 the real 
probability density, Nj = J drrij(r) the number, Aj 
= h 2 (6ir 2 ) 2 / 3 /(2mj), and Vj is the confining trap. Here 
the interspecies coupling is <?i2 = 2-kTi anjmR with the 
reduced mass mn — m\mil(m\ + m^), and an is the 
interspecies scattering length. The interaction between 
intra-species fermions in spin-polarized state is highly 
suppressed due to Pauli blocking term 3Ajn^ /5 and 
has been neglected in Eq. and will be neglected 

throughout. The kinetic energy terms in this equation 
h |V^j'| 2 /(6fJirj') are derived from a hydrodynamic equa- 
tion for the fermions [2(j and contribute little to this 
problem compared to the dominating Pauli blocking term 
in Eq. However, the inclusion of the kinetic energy 
terms in Eq. leads to a smooth solution for the proba- 
bility density everywhere |ll| ■ To keep the algebra simple 
and without losing generality in our calculation we shall 
take equal fermion masses: mi = n%2 = m/3. This sim- 
ulates well the fermion mixtures of ground and excited 
states of 6 Li Q and 40 K Q] atoms observed experimen- 
tally. The Lagrangian density of each fermion component 
in Eq. @J is identical to that used in Refs. [TlL [l3| . 

The mean-field dynamical equations for the system are 
just the usual Euler-Lagrange (EL) equations 
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where Xk, k = 1, 2, 3 are the three space components, and 
j = 1,2 refer to the fermion components. 

With Lagrangian density Q the following EL equa- 
tions of motion are derived in a straight-forward fashion 

MM- 



h 2 vl 
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+ Vj(r) + A 3 n 2 / 3 + g 12 n k % = 0, 
j^k = 1,2. (3) 



This is essentially a time-dependent version of a simi- 
lar time-independent model used recently for fermions 
[Tij . For a stationary state Eqs. © yield the same 
result as the formulation of Ref. [14( . For a system 
with large number of fermions both reduce to lid Il4l 
the well-known Thomas- Fermi approximation |2ll |22| : 
rij = [(jUj — Vj)/Aj] s ' 2 , with fj,j the chemical potential. 

In the spherically-symmetric case we take V}(r) = 
^(3rrij)uj 2 r 2 and ^-(rji) = ipj(r;t). Now transforming 
to dimensionless variables defined by x — y/2r/l, r = tu>, 
I = yJKJJrrujJ) and 
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<t>j(x;r) = 0, 



(5) 



where j ^ k = 1,2 and N n = 3(37r/V,/2) 2 / 3 , Af jk = 
Q^/2Nkav2./L The normalization of the wave-function 
components is given by j Q dx\<j>j(x; t)\ 2 = 1, j = 1, 2. 

We next reduce Eq. to a minimal quasi-one- 

dimensional form in a cigar-shaped geometry where the 
confining trap of anisotropy v has the form Vj(r) = 
^3mjU! 2 (p 2 + v 2 z 2 ), with p the radial vector and z the 
axial vector. For a cigar-shaped geometry v << 1, we 
consider solutions of Eq. (0) of the type 9j(r,t) = 



Njipj(z,t)i/)f\p), where 
\^\p)\ 2 



3rrijUj 



exp 



SmjUjp' 1 
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corresponds to the respective circularly symmetric 
ground-state wave function in the absence of nonlinear 
interactions and satisfies 



2„,,(0) 



6m 



i 2 2 / (0) 

+ 2 m J' w P r) 



fkuip 



(0) 



(7) 



with normalization 2tt J °° \ipj (p)\ 2 pdp = 1. Now the dy- 
namics is carried by <fj(z, t) and the radial dependence 
is frozen in the ground state ip^ (p) . The separation of 
the variables is suggested by the structure of Eq. © . 
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Averaging over the radial mode ipj°\ i.e., multiplying 

Eq. J2J by ipj°^*(p) and integrating over p, we obtain the 
following one-dimensional dynamical equations [TH Il2| : 



in— 



n <9 2 3 2 2 2 
+ —vhav co z 



dt & mj dz 2 ' 2"" J 
F n N 2 / 5 \up^ + F 3k N k \ Vk \ 



i Pj (z,t)=0, (8) 



where j ^ k = 1, 2 and 
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where M12 = "Z>rn\mil (m\ +7712). In Eq. ||SJ| the normal- 
ization is 1^(2, £)| 2 <iz = 1. 

In Eq. JSJ), we introduce the dimensionless variables 
r = tvw/2, x = z/l z , tfij = ^fQfj, with l z — i/S/(rom), 
so that 



^"c^ +c2a;2+A/ ' jfc|</,fc|2 



,4/3 



%( a ;,r)=0, 7^* = 1,2, (11) 



where c = 1, A/}ft = 12ai 2 (Afc/j/)/Z z , and = 
9(67riV : ,/i/) 2 / 3 /5. However, by taking c^l we can mod- 
ify the harmonic trap. In Eq. Ijlll) . the normalization 
condition is given by J_ \(f>j(x,r)\ dx = 1. 




(a) 



-40 -30 -20 -10 10 20 30 40 
z (urn) 




(b) 



-40 -30 -20 -10 10 20 30 40 
z (urn) 



0.4 



sT 0.3 

E 

3 0.2 



JT0.1 




a -IWnm N 1= 1000 

a 12 -l50nm V = 10 



c= 1 




-40 -30 -20 -10 10 20 30 40 
(c) z (urn) 




(d) 



-40 -30 -20 -10 10 20 30 40 



III. NUMERICAL RESULTS 

We solve the coupled mean-ficld-hydrodynamic equa- 
tions JSJ and Ijlll) numerically using a time-iteration 
method based on the Crank-Nicholson discretization 
scheme elaborated in [23j. We discretize the mean-field 
equations using time step 0.001 and space step 0.05. In 
the quasi-one-dimensional case we span x from —60 to 
60, and in the three-dimensional case we span x from 
to 30. In both cases we start with the Gaussian 
ground-state solution of the linear harmonic oscillator 
problem as input at t — while we set the nonlincaritics 
JVjk = Afjj = 0. In the three-dimensional case governed 
by Eq. we take r = 0) = (tt/2)- 1 / 4 exp(-a; 2 /4) 
and in the one-dimensional case governed by Eq. I|ll|) we 
take 4>j(x, t = 0) = 7r -1 / 4 cxp(— x 2 /2). With these initial 
states we perform time evaluation of Eqs. JSJ and (|llfl . 
During the course of time evolution the nonlinear terms 
Afjk and Afjj are slowly introduced at the rate of 0.0001 
at each time step of 0.001. The solution evolves slowly 
with time and attains the final stationary configuration 



FIG. 1: (Color online) The fermionic profiles <fij(z) of Eq. JSJ 
vs z for N\ — 1000 and N2 = 10 in the quasi-one-dimensional 
case {v = 0.1) for ai2 = (a) 50 nm, (b) 100 nm, (c) 150 nm, 
and (d) 400 nm. 

after the full desired nonlinearity is introduced slowly. 
It is important to introduce the nonlinearity slowly to 
achieve a final stationary configuration. 

In our numerical investigation we take / = 1 pm and 
consider the fermions to have the mass of 40 K atoms 
and to occupy two levels of this fermionic atom. In the 
spherically-symmetric case this corresponds to a radial 
frequency oj « 2tt x 83 Hz. Consequently, the unit of 
time is 1/u w 2 ms, and the unit of length l/y/2 ~ 0.7 
pm. In the quasi-one-dimensional case we take v = 0.1 
and l z = 1 /im corresponding to an axial frequency 
ljv Ri 27r x 83 Hz. Consequently, the unit of time is 
2/(yu>) w 4 ms, and the unit of length l z = 1 pm. 

After a small experimentation it is realized that both 
in Eqs. JSJ and (|11|) . demixing is most favored when N\ 
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FIG. 2: (Color online) The fermionic profiles <j>j{x) of Eq. © 
vs. x in dimensionless units for N\ = 10000 and N% = 100 in 
the spherically-symmetric case for ai2 = (a) 50 nm, (b) 300 
nm, (c) 400 nm, and (d) 800 nm. 



and N2 are widely different from each other. Under this 
condition one of the components is much larger in size 
than the other and easily expels the smaller component 
from the center of the trap upon an increase in inter- 
species repulsion or trap frequency and we shall consider 
only this asymmetric configuration here. In previous in- 
vestigations |l8j| demixing in the symmetric configuration 
(Ni = N 2 ) was considered. The interspecies repulsion 
can be manipulated through the scattering length a\2 
by varying a background magnetic field near a fermion- 
fermion Feshbach resonance [24j . 

First, we present results for the quasi-one-dimensional 
case from a solution of Eq. Ijlll) . After some experimen- 
tation we take in the initial DFFM N x = 1000, N 2 = 10, 
and a\2 = 50 nm, v = 0.1. This corresponds to nonlinear- 
ities Mi » 5917, M2 ~ 274, M2 = 60 and Mi = 6000. 
The two overlapping Gaussian-type fermionic profiles in 
this case are shown in Fig. 1 (a). Next we increase the 
interspecies repulsion by increasing the interspecies scat- 



tering length ayi- With the increase of repulsion the com- 
ponent with smaller number of fermions is slowly deiven 
out of the central region which continues to be solely oc- 
cupied by the component with larger number of fermions. 
In Figs. 1 (b), (c) and (d) we plot the results for ayz — 
100 nm, 150 nm, and 400 nm, from where we find that 
a complete demixing has occured for 012 = 400 nm. In 
the last case two symmetrical peaks of the component 
with smaller number of fermions appear where the den- 
sity of the component with larger number of fermions is 
practially zero. The interesting thing to note is that in 
Figs. 1 (a) — (d), the profile of the degenerate gas with 
larger number of fermions remains practially unchanged, 
while the weaker component with much smaller number 
of fermions is expelled from the central region with the in- 
crease of repulsion. The mutual repulsion has practically 
no effect on the "heavier" component when it expells the 
"lighter" component from the central region. Had we 
considered the symmetric case with N\ = N2 = 1000 a 
much larger a\2 would be necessary for a complete demix- 
ing. 

Next we consider demixing in the three-dimensional 
spherically-symmetric case from the solutions of Eq. JSJ. 
In this case we take in the initial DFFM N x = 10000, 
N2 = 100, and a\2 = 50 nm. This corresponds to non- 
linearities Mi ~ 3914, M12 ~ 42, Mi ~ 4243, and 
M2 ~ 181. The overlapping fermionic quasi-Gaussian 
profiles of the DFGs in this case are shown in Fig. 2 (a). 
Now we increase the interspecies repulsion by increasing 
the scattering length ai2 • With the increase of repulsion 
among the two species, smooth Gaussian profiles of the 
DFFM are destroyed and the component with smaller 
number of fermions is slowly driven out of the central 
region which continues to be occupied by the component 
with larger number of fermions. In Figs. 2 (b), (c), and 
(d) we plot the results for 012 = 300 nm, 400 nm, and 800 
nm, respectively. From Figs. 2 we see that demixing has 
increased with the increase of repulsion among the two 
components. In Fig. 2 (d) the component with smaller 
number of fermions is practically located outside the re- 
gion occupied by the component with larger number of 
fermions. 

The dimensionless parameter controlling the validity 
of the present dilute (weakly interacting) mean-ficld- 
hydrodynamic model is the number of particles in a scat- 
tering volume I a, 1 3 , where a is the scattering length. This 
can be written as n|a| 3 , where n is the average density 
of the gas |2l|. When n\a\ 3 « 1 the system is said 
to be dilute or weakly interacting. The opposite limit 
n|a| 3 >> 1 is known as the unitarity limit. It is interest- 
ing to find out if demixing in the three-dimensional case 
shown in Fig. 2 occured in the dilute limit. In Fig. 2 (c) 
demixing is visible for \a\ = 400 nm where there are 10000 
atoms in a sphere of radius R » 18 /xm. Consequently, 
n = 10000/(47ri? 3 /3) w 0.4 ^m" 3 and \a\ 3 w 0.064 /im 3 , 
so that n\a\ 3 « 0.025. So demixing starts in the dilute 
weakly interacting limit and continues in the unitarity 
limit. 
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FIG. 3: (Color online) The fermionic profiles <Pj{z) of Eq. JSJ 
vs. z for Ni = 1000, N 2 = 10 and a 12 = 100 nm in the 
quasi-one-dimensional case (v = 0.1) for c = (a) 3 and (b) 8. 



The present condition of demixing is consistent with 
that of Ref. , where they have reported demixing as 
a function of the parameter 22] 



pot 



(12) 



with a = 2 1 / 3 3 1 / 6 (8192/2835tt 2 ) « 0.44 and AT the num- 
ber of atoms, -Ej n t is the interaction energy and -Ep t 
is the potential energy of the harmonic oscillator trap. 
The parameter 7 compares the interaction energy with 
the potential energy in the harmonic trap and provides a 
more reliable account of the words "weakly interacting" . 
A small value of the parameter 7 together with the condi- 
tion n\a\ 3 << 1 is necessary for the DFFM to be weakly 
interacting. In the situation of Fig. 2 (c) (a = 400 nm) 
demixing has started for 7 w 0.44 x 10000 1/6 x 0.4 w 0.8 
and is absent in the situation of Fig. 2 fb_) (a = 300 
nm) where 7 ~ 0.6. In the study of Ref. yj| demixing 
started at 7 ss 0.7, and was absent for 7 w 0.54 quite con- 
sistent with present finding. The small value of t(< 1) 
obtained in the present study as well as in Ref. [lg con- 
firm the dilute or weakly interacting nature of the DFFM 
in the demixing limit. However, according to the condi- 
tion obtained in Ref. |l7j . demixing should occur only 
in the unitarity limit. We do not know the reason for 
this discrepancy. Further investigations are needed for 
clarification. 

Now we investigate the possibility of demixing in 
a DFFM by increasing the strength of the harmonic 
trap. We illustrate the results only in the quasi-one- 
dimcnsional case described by Eq. (jTT)l . From Fig. 1 



(b) we find that there is no real demixing in the solution 
of Eq. JTTJ for N x = 1000, N 2 = 10, a 12 = 100 nm and 
c = 1. In this case we find that demixing appears as we 
increase c keeping other parameters fixed. The results 
are illustrated in Figs. 3 (a) and (b) for c = 3 and 8. 
The increase of parameter c increases the trap frequency 
in the same proportion. From Fig. 3 (a) we find that 
demixing already appears for c = 3 and increases in Fig. 
3 (b) for c = 8. For a large c, the appearance of demixing 
for a smaller a\i can be understood from the structure 
of Eq. Hll|) . An increase of the trapping frequency cor- 
responds to a decrease in the harmonic oscillator length 
l z . Then a smaller value of 012 can lead to a larger non- 
linearity Mjk as in the definition of A/}fe, a\2 is always 
scaled with respect to l z . However, with the increase of 
the trapping frequency the scaling length l z is reduced 
and the degenerate gas is confined to a smaller region in 
space as can be seen from Figs. 3 (a) and (b). 



IV. SUMMARY 

Using a coupled set of time-dependent mean-field- 
hydrodynamic equations for a trapped DFFM derived 
using the Lagrangian density for a DFG suggested re- 
cently |lll Il2j |. we study the phenomenon of demixing in 
a trapped DFFM. We show that, for a repulsive inter- 
species fermion-fermion interaction, demixing is favored 
in a DFFM with different number of fermions. Demix- 
ing appears for a large value of interspecies scattering 
length ayi denoting a large repulsion between the com- 
ponents. As the interspecies repulsion increases, the 
fermionic component with smaller number of fermions 
is expelled from the central region which is solely occu- 
pied by the component with larger number of fermions. 
Demixing in a trapped DFFM can be realized experimen- 
tally by changing the interspecies scattering length ayi 
to a large positive value by exploiting a fermion-fermion 
Feshbach resonance [24] . Demixing also appears with the 
increase of trapping frequency for a smaller a\2- This 
later possibility to achieve demixing is more attractive 
experimentally as it is easier to increase the trapping fre- 
quency than to increase a\2 near a Feshbach resonance 
by manipulating a background magnetic field |24|. A 
proper treatment of DFFM should be performed using 
a fully antisymmetrized many-body Slater determinant 
wave function [15j as in the case of atomic and molecu- 
lar scattering involving many electrons |25|. However, in 
view of the success of the hydrodynamic model in other 
contexts 0, 0, we do not believe that the present 
study on mixing-demixing in a DFFM to be so peculiar 
as to have no general validity. 
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